Family of exactly solvable models with an ultimative quantum paramagnetic ground state 
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We present a family of two-dimensional frustrated quantum magnets solely based on pure nearest-neighbor 
Heisenberg interactions which can be solved quasi-exactly. All lattices are constructed in terms of frustrated 
quantum cages containing a chiral degree of freedom protected by frustration. The ground states of these models 
are dubbed ultimate quantum paramagnets and exhibit an extensive entropy at zero temperature. We discuss the 
unusual and extensively degenerate excitations in such phases. Implications for thermodynamic properties as 
well as for decoherence free quantum computation are discussed. 
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The search for exotic phases of strongly correlated quantum 
matter possessing unusual physical properties is a fascinat- 
ing and active research area. Exactly solvable models played 
an essential role toward this end. Important examples like 
the Heisenberg chain or the frustrated Majumdar-Ghosh point 
with free or gapped spinon excitations exist in one dimension 
[1-3]. In two dimensions, the frustrated Shastry-Sutherland 
model constitutes the paradigm example for a valence bond 
solid [4]. Its experimental realization in SrCu2(BC>3)2 has 
spurred intense activity. Finally, the exactly solvable toric 
code [5, 6] and Kitaev's honeycomb model [7] have become 
standard models rigorously exhibiting topologically ordered 
phases. These have received enormous attention in the con- 
text of topological quantum computation. 

As already indicated above, one driving knob to tune quan- 
tum fluctuations is geometrical frustration. Indeed, the most 
promising candidates for experimental realizations of a quan- 
tum spin liquid (SL) are strongly frustrated systems like the 
so-called Herbertsmithite [8, 9] or the undoped parent com- 
pound k-(BEDT-TTF)2Cu2(CN)3 of the organic superonduc- 
tors [10]. In Herbertsmithite, the geometrical frustration in the 
Heisenberg model on the kagome lattice is proposed to stabi- 
lize a SL state [11], while multi-spin interactions are relevant 
for the organic compound [12]. 

Theoretically, different quantum SL possessing very differ- 
ent physical properties have been proposed. Gapped SL have 
been exactly demonstrated in the quantum dimer model on 
non-bipartite lattices [13, 14] and in the toric code [5, 6]. Fur- 
ther interesting examples are the U(l) critical SL discussed 
in the framework of the kagome model [11] and Bose-metal 
phases discussed for triangular topologies with multi-spin in- 
teractions [15]. Finally, there are intriguing proposals for 
so-called ultimate co-operative paramagnets in context of the 
transverse field Ising model on the kagome lattice [16]. In all 
these cases, no rigorous examples of such exotic phases in two 
dimensions are known. Here, we present a family of quasi- 
exactly solvable microscopic models having quantum disor- 
dered, extensively degenerate ground states and ultra short- 
ranged spin correlations. We dub these distinct class of SL 
ultimate quantum paramagnets (UQP). 

Our models might well be realized in experiment and, addi- 
tionally, could be relevant for quantum information (QI). First, 
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FIG. 1: (Color online) (a) Heisenberg model with exchange J on 
a periodically coupled chain of Nt triangles, (b) Illustration of the 
two degenerate ground states. Thick (green) bonds represent sin- 
glet states, (c) The shortest segment of a Shastry-Sutherland chain 
consists of two vertical and one horizontal dimer. Note that the two 
dimers on the edges are part of the quantum cage. 



the models consist solely of pure nearest-neighbor spin- 1/2 
Heisenberg antiferromagnets (HAF), i.e. coupled sites i and 
j interact via Si ■ Sj. This kind of exchange is realized in 
many physical systems. In the QI context, we are especially 
motivated by coupled quantum-dot systems [17-19], possible 
patterning studies of quantum cellular automata (QCA) [20], 
trapped ions [21, 22], and by using triangular magnetic clus- 
ters like CU3 [23, 24]. Suitably designed Josephson junction 
arrays (JJA) are also of interest here, since the JJA can be 
mapped onto S = 1/2 Heisenberg-like models [25]. We draw 
attention to the fact that the core element in QCA computation 
is a bistable "cell" capable of interacting with its neighbors, 
which is precisely the strategy used here. Indeed, protected 
chiral quantum bits as also engineered in our case in a scalable 
fashion have been the focus of recent works in QI [23, 24]. 

Let us first introduce the construction recipe for the lattices 
we study in this work. This is crucial since all interesting 
physical properties emerge from the special frustrated topol- 
ogy [26] . The elementary building block is a finite chain of N t 
coupled triangles with periodic boundary conditions, as illus- 
trated in Fig. la. Clearly, the ground state of this model is two- 
fold degenerate for Nt > 2. These states are singlet coverings 
which are either oriented to the right or the left (see Fig. lb). 
Thus, this two-level system represents a chiral pseudo-spin. 

Next, we put this chain segment on a ring as shown in Fig. 2 
for Nt € {3, 4, 6}. Now we add 7Y t dimers on the boundary 



to form cage-like structures. Since the boundary dimers are 
coupled to the inner triangles in a Shastry-Sutherland fashion, 
singlets are formed on the outer dimers in the ground state 
of the full cage and the two-fold chiral degeneracy remains 
protected. Hence, we dub this unit a chiral quantum cage. 

The final step is to couple such cages. To this end, we in- 
troduce one additional dimer between two cages such that the 
inter-cage coupling involves the shortest segment of a one- 
dimensional Shastry-Sutherland chain (see Fig. lc) as illus- 
trated in the right side of Fig. 2. The boundary dimers are 
coupled to the connector with strength J' = x J. In this way, 
one can build any two-dimensional lattice of coupled cages. 
We are only aware of two works dealing with Hubbard [27] 
and Heisenberg [28] models on decorated lattices in the con- 
text of unconventional order, but gapped spin liquid ground 
states have never hitherto been considered in this context. 

We want to study the ground states and elementary excita- 
tions as a function of x — J' J J. Our focus will be the regime 
x < x c where, remarkably, the ground states of these models 
can be found exactly. To see this, notice that the full lattice can 
be covered by singlet configurations on each cage. The inter- 
cage coupling involves the basic Shastry-Sutherland unit, i.e. 
two vertical singlets from neighboring cages are linked by one 
horizontal singlet. Since, as a result, the chiral degree of free- 
dom on neighboring cages are effectively decoupled, the exact 
ground states for x < x c are a direct product of the dimer cov- 
erings on each of the N c cages. In contrast to the dimer solid 
in the Shastry-Sutherland model [4], the crucial point here is 
that this number is extensively large. The ground state degen- 
eracy scales like 2 Nc , leading to a finite entropy ln(2) per cage 
at T = in the thermodynamic limit. As x increases, a phase 
transition at x = x c will certainly occur (see below). 

But let us first discuss the elementary excitations in the 
UQP ground state for x < x c . They can also be determined 
quasi-exactly. One finds three sorts of massive particles with 
total spin one above the infinitely degenerate ground state: i) 
excitations on the connector, ii) excitations on the boundary 
dimers, and iii) excitations inside the cage (stars). 

Excitations on the connector are localized triplets. This is 
a direct consequence of the peculiar Shastry-Sutherland ge- 
ometry in one dimension which surpresses any fluctuations to 
neighboring dimers. The connecting dimer can be only in a 
singlet or a triplet state, both eigenstates of the full problem. 
The total Hilbert space separates into different blocks belong- 
ing to a different set of singlet-triplet configurations. A single 
triplet costs the exact energy J independent of x. Clearly, the 
true low-energy excitations of the model belong to the sector 
with singlets on all connecting dimers. 

The second sort of excitations are on dimers on the bound- 
ary of a cage. These excitations also stay local. But triplets 
are not exact eigenstates, since they can virtually create and 
annihilate triplets on the inner star as well as on the attached 
connector, thereby reducing its energy. These excitations are 
local triplons [29] and constitute the excitations with lowest 
energy. This mode is displayed in Fig. 3 for N t = 3. One finds 
a gap 0.23J at x — 0. The gap is lowered as x increases and 
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FIG. 2: (Color online) Three representative lattice models with (a) 
N t = 3, (b) N t = 4, and (c) N t = 6. The left sides illustrates the 
chiral quantum cage of each lattice which contains a triangle chain 
shielded by surrounding dimers. The right sides represent 2d lattices 
made by adding dimers which couple cages. The dimer connecting 
two cages forms the shortest segment of a Shastry-Sutherland chain 
(illustrated in black). 

the mode crosses the low-energy chiral manifold at x « 0.7. 
This signals a phase transition as long as no other levels cross 
before. Below, we argue this to be the case here. 

Finally, one has the excitations inside the star of a cage. The 
lowest modes correspond to very weakly dispersive, massive 
triplons. Let us note that these excitations do not depend on x 
and can therefore be identified as the horizontal gaps in Fig. 3 
for Nt = 3. Studies on the infinite chain Nt = oo found 
completely local excitations with an energy gap 0.22J [30]. 
For finite Nt the energy gap is slightly larger. One finds 0.32 J 
for N t = 3. The weakly dispersing character is nicely under- 
stood from the dimerized limit of the inner star of a cage by 
introducing a modulation J ± S on the spikes of the triangles 
(see Fig. 1). In the limit S = 1 the ground state is the product 
state of singlets and excitations are local triplets. Turning on 
a finite 6, a first hopping processes of the triplet occurs in per- 
turbation order J 2 ^' -1 ) . One first excites all triplets of the 
star one after the other (Nt — 1 operations) and subsequently 
annihilates Nt — 1 triplets such that the triplet has effectively 
hopped. Thus, one directly understands that the excitations 
become increasingly local with increasing 7Y t , 
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FIG. 3: (Color online) Localized excitations on a chiral quantum 
cage with N t — 3 as a function of x. Thick grey line at energy J 
corresponds to the exact triplet on the connecting dimer. Blue hori- 
zontal lines represent excitations inside the star. Dashed lines involve 
localized triplons on the outer shell which are the only single modes 
depending on x. The lowest of these modes is the energy of a single 
triplon. Other dashed lines corresond to a triplon plus an excitation 
inside the star. Number in square brackets gives the degeneracy on 
the cluster shown in the inset which has been diagonalized. 

A fourth exactly known class of excitations is the combined 
creation of a triplon on the outer shell plus excitations of the 
star. These modes certainly depend on x due to the presence of 
the local triplon. Interestingly, a strong attraction is found for 
certain modes which leads to the presence of levels at rather 
low energy (see Fig. 3) [37]. 

In sum, we have seen that ground states and all relevant 
low-energy excitations for x < x c can be found quasi-exactly. 
Excitations on the outer shell of the cages have the lowest 
energy. We stress that these excitations are localized on the 
quantum cage. This implies an extensive degeneracy for these 
states as well. Since ground state correlations are strictly zero 
for distances larger than the extension of a quantum cage, we 
have found a family of exactly solvable models with UQP 
ground states. 

This can actually be extended in two respects. First, one can 
change the stars to different geometries also possessing a two- 
fold degenerate ground state, e.g. a fully-frustrated plaquette 
with the minimal number of four spins. If one couples the 
fully frustrated plaquettes in a completely analogous fashion 
as before, one again finds exactly an UQP ground state for 
x < x c (see Fig. 4a). 

A second route is to change the nature of the connector. 
Indeed, a direct coupling x J of stars by vertical dimers also 
protects the UQP ground state for x not too large. A graphical 
illustration for the case JV t = 3 is shown in Fig. 4b. Again, 
excitations inside the star and the local triplon centered on the 
connecting dimer will be exact excitations of the full problem. 

As x increases, a phase transition will certainly occur. De- 
tails depend on the nature of the connector. For the case of a 
direct star coupling just discussed this will likely be the triplon 
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FIG. 4: (Color online) a) The fully-frustrated plaquette is the smallest 
inner star, b) A direct coupling of stars (Nt = 3 shown) leads to a 
different family of models, c) Illustration of the two possible close 
packings of JVt =6 quantum cages on the Kagome lattice. Left 
figure represents the proposed valence bond solid ground state with 
the 36-site unit cell. Right figure represent a diluted Kagome lattice 
where sites marked by grey squares are removed. 

mode on the connector. If true, a magnetically ordered triplon 
condensate will be stabilized for large x. This triplon mode 
can be exactly captured by studying a cluster of two neighbor- 
ing stars coupled by a single dimer. We find a critical value 
x c ps 0.75 for N t = 3. 

In contrast, for the connector comprised of the Id segment 
of a Shastry-Sutherland chain discussed first, a different sce- 
nario is likely. One knows for the Shastry-Sutherland model 
that triplons strongly attract in the sector with 5 = when 
placed on neighboring vertical dimers, meaning the two ends 
of a connector in our case [31]. We find strong evidence for 
this scenario for coupled fully-frustrated plaquettes. Here, one 
observes on a cluster of two plaquettes and one connector that 
the UQP becomes unstable for x c ps 0.85. Instead of a single 
triplon, it is again a singlet state keeping the ground state de- 
generacy of the UQP phase. The two triplons on neighboring 
cages bind into a singlet bound state. 

We stress that the above considerations do not fix the nature 
of the phase transition. It might well be that the transition for 
both cases turn out to be first order as it is known for the one- 
dimensional Shastry-Sutherland chain [32]. If so, the phase 
transition point would tend slightly to lower values of x. 

Next, we want to discuss a connection of the case N t = 6 
with Shastry-Sutherland connector to the strongly debated 
Heisenberg model on the Kagome lattice. There are two pos- 
sible close packings of chiral quantum cages on the Kagome 
topology (see Fig. 4c). The first packing is such that honey- 
combs remain between the cages. Assuming a spontaneous 
dimerization, one recovers the valence bond solid consistent 
with 36-site unit cell proposed to be linked to the true ground 
state of the Heisenberg model on the Kagome lattice [33, 34]. 
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The second packing is denser but leaves dangling spins which 
is clearly unpreferable. But removing these spins yields a 
depleted Kagome lattice. Quite remarkably, one can show 
that the depleted model has VBS order with a 24-site unit 
cell. The easiest way to see this is by identifying the chi- 
ral pseuodspin states as exact eigenstates of a single cage. In 
2nd order inter-cage degenerate perturbation theory, the effec- 
tive pseudospin model turns out to be precisely an Ising fer- 
romagnet in a transverse "magnetic" field (see also Ref. 35): 
Heff = -J T~T* - h T£, with j,h « J . The 
"field" stabilizes a VBS state with "ferro" chiral order. The 
elementary excitations are therefore weakly dispersive "pseu- 
dospin" flips involving co-operative flipping of the six dimers 
on the internal star. This succinctly illustrates the crucial in- 
fluence of the unusual ground states found here on geometry, 
with disordered valence bond liquid, valence bond solid and 
magnetically ordered phases arising for differing geometries. 

Let us finally discuss the specific heat. For x < x c , 
the decoupled (disordered) phase is characterized by a spe- 
cific heat entirely due to the fluctuating chiral doublet, giv- 
ing C(T) ~ T" 2 , reflecting the finite T=0 entropy per 
site. In the case where magnetic order arises for x > 
x c , quenching of this entropy will yield C(T) ~ T 2 . For 
the depleted Kagome magnet, two weakly dispersive bands 
of singlets, defined by = ±(V e fc + ~ h ) with 

efc = 2j [cos k y + 2 cos (y/3k x /2) cos (k y /2)] and = 
(2/3)jsm(k y /2)[cos(V3k x /2)~cos(ky/2)] for h « j, 
give rise to phonon-like excitations with velocity v ~ j, and, 
in two dimensions, this will yield C(T) ~ exp(—h/t) for 
T -> and C(T) ~ T 2 for higher temperatures. 

In conclusion, we have presented families of frustrated 
quantum magnets with hitherto undiscovered gapped SL 
ground states. The ground state as well as elementary exci- 
tations can be determined exactly and both possess a macro- 
scopic degeneracy at T = 0. One important question is the 
fate of the finite entropy upon heating these systems. Indeed, 
thermal fluctuations induce a finite triplon density on each 
cage. These triplons correspond to different effective pseudo- 
spins (not to be confused with the chiral pseudo-spins T M ). 
This leads to an inter-cage coupling, producing an ordered 
phase continuously connected to that found at T = for x 
> x c . 

For experimental realizations, we only demand flexibility 
on the lattice design since our models solely contain nearest- 
neighbor Heisenberg couplings. We suggest that experimental 
realizations of QCA in suitable geometries, artificially fabri- 
cated JJA, or trapped ions may exhibit some of the unusual 
features found in this work. Arrays of magnetic clusters, e.g, 
of CU3 triangles suitably engineered on a substrate, may also 
reveal parts of the exotica proposed here. In particular, given 
the protected chirality, an exciting option for future study is to 
use this aspect in a setting favorable for quantum computation 
in "interaction free subspaces" [36] to minimize the trouble- 
some decoherence problem in quantum information. 
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